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We calculate the next-to-leading-order (NLO) QCD corrections to e + e~ — > J /ipgg via color singlet 
J/tp{ 3 Si) at the B factories. The result shows that the cross section is enhanced to 0.373 pb by a 
K factor (NLO/LO) of about 1.21. By considering its dependence on the charm quark mass and 
renormalization scale, the NLO cross section can range from 0.294 to 0.409 pb. Further including 
the if) 1 feed-down, a(e + e~ — » J/ipX(non cc)) is enhanced by another factor of about 1.29 and 
reaches 0.482 pb. In addition, the momentum distributions of J/ip production and polarization are 
presented. Recent measurements from Belle agree well with our prediction for the cross section and 
momentum distribution. It is expected that this process can serve as a very good channel to clarify 
the J/ip polarization puzzle by performing further experimental measurements. 
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Perturbative quantum chromodynamics are successful 
to describe large momentum transfer processes at quark 
level due to its asymptotic freedom property. But it falls 
into the nonperturbative range for quark hadronization 
to form the final state hadrons which are measured in 
experiments. Therefore the quark hadronization is usu- 
ally described by phenomenological models and quite far 
away from the first principle QCD Theory. However, in 
the heavy quarkonium case, a naive perturbative QCD 
and nonrelativistic treatment for the bound state is ap- 
plied straightforwardly to the related decay or produc- 
tion processes. It is called color-singlet mechanism. To 
describe the huge discrepancy of the high-p t J/ip produc- 
tion between the theoretical calculation based on color- 
singlet mechanism and the experimental measurement 
at Tevatron, a color-octet mechanism I] was proposed 
based on the non-relativistic QCD (NRQCD) 2]. It al- 
lows consistent theoretical predictions to be made and 
to be improved systematically in the QCD coupling con- 
stant a s and the heavy-quark relative velocity v. In re- 
cent years, there is a huge data collection in the B fac- 
tory experiments. Based on that, many J/ip production 
processes were observed [1, 0, [j| in the past. Now the 
integrated luminosity is more than 850 fb _1 at the Belle 
detector at the KEKB and it is about 20 times larger than 
the integrated luminosity 32.4 fb _1 , based on which the 
inclusive J / ip production was measured [3J, |J] . Therefore 
it supplies a very important chance to perform system- 
atical study on J/ip production both theoretically and 
experimentally. 

The measurements for exclusive J/ip productions 
e + e _ — > J/ipT) c , J/ipJ/ip, J/ipXcj at the B factories 
have shown that there are large discrepancies between 
the leading-order (LO) theoretical predictions 
in NRQCD and the experimental measurements [4|,| 
It seems that such discrepancies can be resolved by intro- 
ducing higher order corrections @, [ll|, [HI, [l3|, [l4[ : next- 
to-leading-order (NLO) QCD corrections and relativistic 
corrections. 



The cross section for inclusive J/ip production in e + e 
annihilation was measured by BABAR [l5| as 2.54 ± 
0.21 ± 0.21 pb and Belle @, [| as 1.45 ± 0.10 ± 0.13 pb. 
These measurements include both J/ip + cc + X and 
J/ip + X(non cc) parts in the final states. Many the- 
oretical studies @, [13, OS OS HI HH IH [H, [Hf have 
been performed on this production at LO in NRQCD 
and the results for inclusive J/ip production cover the 
range 0.6 ~ 1.7 pb depending on parameter choices. A 
further analysis by Belle |4j gives 
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J/ipcc + X) = 0.87±g;?g± 0.17 pb. (1) 



It is about 5 times larger than the LO NRQCD predic- 
tion However, this large discrepancy was partially 
resolved by considering both NLO correction and feed- 
down from higher excited states [l6| ■ The above measure- 
ments infer that a[e + e~ — > J/ip + A"(non cc)] = 0.6 pb. 
For this part, the contributions from the color-singlet 
and color-octet contributions for the processes, e + e~ — > 
J/V (1) ( 3 5i).gg, J/7/' (8) ( 1 5o, 3 P/)g, are about 0.2 pb and 
0.27 pb, respectively, at the LO in NRQCD 24]. How- 
ever, the signal of the color octet was not found in the 
experiment 0, flBI ]. Therefore, the experimental mea- 
surement by Belle is about 3 times larger than the theo- 
retical prediction from color singlet at LO, and can be 
much more than 3 time s by BABAR. Recently many 

studies [II OS OS OS [Mill III haVC shown that hi g her 
order corrections within NRQCD framework are very im- 
portant. To achieve a reasonable theoretical prediction 
for the process e + e~ — > J/ipgg, in this letter, we present 
a NLO QCD calculation to this process. 

The related Feynman diagrams which contribute to 
the LO amplitude of the process e + (pi) + e~(p2) — ► 
J/tp(P3) +<?(P4) +5(^5) are shown in Fig. [TJa), while the 
others can be obtained by permuting the places of the 
virtual photon and gluons. In the nonrelativistic limit, 
the NRQCD factorization formalism is used to obtain the 
total cross section in n — 4 — 2e dimensions as 
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FIG. 1: Typical Feynman diagrams, (a): LO; (b) e + e — * J/ipgqq; (c) e + e 
the counter-term diagrams, including corresponding loop diagrams. 



J/tpggg; (d): One-loop. Groups (di) — (ds) are 
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by introducing a dimensionless kinematic variable s = 
s/(2m c ) 2 , where s is the squared center-of-mass energy, 
e c and m c are the electric charge and mass of the charm 
quark, respectively, and /3 = y/ s(s — 1). R s (0) is the 
radial wave function at the origin of J ji>. The approxi- 



mation M 



2m c is taken. 



At NLO in a s , there are virtual corrections which arise 
from loop diagrams. Dimensional regularization has been 
adopted for isolating the ultraviolet (UV) and infrared 
(IR) singularities. UV-divergences from self-energy and 
triangle diagrams are canceled upon the renormalization 
of QCD. Here we adopt the same renormalization scheme 
as Ref. 12] . There are 111 NLO diagrams in total, includ- 
ing counter-term diagrams. They are shown in Fig. [TJc?) , 
and divided into 8 groups. Diagrams of group (d$) con- 
tain Coulomb singularities, which can be isolated and 
mapped into the cc wave function. Although the Feyn- 
man diagrams are similar, the calculation of tensor and 
scalar integrals is much more complicated than that in 
Ref. [25j], because there is one more variable, which is the 
mass of the virtual photon, in this calculation. Again, 
the calculation was done automatically with our Feyn- 
man Diagram Calculation package (FDC) 26]. 

The real corrections arise from two processes, e + e~ — > 
J/ipgqq and e + e~ — * J/ipggg. The related Feynman dia- 
grams for these two processes are shown in Fig. [TJ6) and 
0Jc). The phase space integration for them will generate 
IR singularities, which are either soft or collinear and can 
be conveniently isolated by slicing the phase space into 
different regions. We adopt the two-cutoff phase space 
slicing method [29] to decompose the phase space into 
three parts by introducing two small cutoffs, 5 S and 5 C . 
And then the real cross section can be written as 

a = a + a +cr, \6) 

where a s from the soft regions contains soft singularities 



and is calculated analytically under soft approximation. 
It is easy to find that the soft singularities for a gluon 
emitted from the charm quark pair in the S-wave color 
singlet J/ip are canceled by each other. And we have 
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where fj, is the renormalization scale and 6 g is the angle 
between two gluons in the pi + p 2 rest frame. a H from 
the hard collinear regions contains collinear singularities 
and can also be factorized. Here we have 



da HC = da^ 
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where 5s = S s /[1 — (ps +Pj) 2 /s] and riif is the number 
of active light quark flavors. The hard noncollinear part 
a HC is IR finite. Finally, all the IR singularities are 
canceled analytically. After adding all the contribution 
together, the cross section at NLO can be expressed as 
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where (3q is the one-loop coefficient of the QCD beta func- 
tion. 
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To study the polarization of J/ip production, we define 
the angular distribution A and polarization factor a as: 



d 2 a 



d cos 9dp 



S(p)[l +A(p) cost 



a 



<7j> 



(7) 



where p and 9 are the 3-momentum and production an- 
gle of J/ip in the laboratory frame, gt and <tl are the 
transverse and longitudinal polarized cross section. To 
calculate a, we use the same method to represent the 
polarized cross section as Eqs. (8) and (9) in Ref. [25|. 
This method is found numerically unstable in a small re- 
gion of phase space due to the cancellation of large num- 
bers. Therefore, the momentum distributions for A and 
a contain potentially large numerical errors in our calcu- 
lation for p < 0.5 GeV or p > 4.2 GeV. As regards the 
total cross section and momentum distribution of J/ip 
production, a simplified method is used to calculate the 
amplitude square with very good behavior in numerical 
calculations. But it cannot be applied to the calculation 
of A and a. 

The values of a s and the wave function at the origin 
of J/ip_in the NLO calculation are taken the same as in 
Ref. [13| . The numerical results are showed in Table. [U 
The scale dependence of the cross section is shown at 



m c (GeV) 




<7^(pb) 


a(s) 


<j (1) (pb) 




1.4 


0.267 


0.341 


2.35 


0.409 


1.20 


1.5 


0.259 


0.308 


2.57 


0.373 


1.21 


1.6 


0.252 


0.279 


2.89 


0.344 


1.23 



TABLE I: Cross sections with different charm quark mass m c 
where the renormalization scale fi — 2m c and y^i = 10.6 GeV. 



Fig. [2] and it does improve significantly at NLO. The 
final numerical result can be expressed as 



o- (1) = 0.373 



+0.036 
-0.079 



pb 



(8) 



where the theoretical uncertainty is from the choices of 
m c and [A, with m c = 1.4 GeV and fi — 2m c for the 
upper boundary and m c — 1.6 GeV and /i = y/s/2 
for the lower boundary. The momentum distribution of 
J/ip production are shown in Fig. [3j To included the 
ip' contribution in the momentum distribution with a 
suitable kinematic treatment, m c — m^//2, /i = m*, 
Br(>' -► J/ip + X) = 0.574 and T{ip' -► ee) = 2.19 KeV 
are used. In Fig. [4] the momentum distributions of the 
polarization factor a and the angular distribution coef- 
ficient A of J/ip are shown. Both a and A have slight 
changes at NLO. Furthermore, if the contribution from 
ip 1 is included, the two curves change very little. In ad- 
dition, we find that da/dcos9 is a constant within quite 
large numerical error. 

To further discuss the renormalization scale \x depen- 
dence, two ways are applied with m c = 1.5GeV. In 
one way, to find the maximum cross section by using 
dcr^/dji — 0, <t2L = 0.41 pb is obtained with fx = 
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FIG. 2: Cross sections as function of the renormalization scale 
fj, and the center-of-mass energy of e + e~ yOs. 
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FIG. 3: Momentum distribution of J/ip production with m c 
1.5 GeV and /i = 2m c . 



1.7GeV, which means the calculation is right in the per- 
turbative range and gives an uncertainty of about 0.04 pb 
for the default estimate. In another way, we apply the 
Brodsky, Lapage, Backenzie (BLM) scale setting[3(| and 
find the unique scale choice fi* = 0.74GeV with the cross 
section crW(/z*) = a^(fi*)(l - 9.1a s {fi*)/n). It means 
that BLM way can not give good convergence behavior 
in this case. 

In summary, we have calculated the NLO QCD correc- 
tion to e + e~ — > J/ipgg at the B factories. It increases the 
cross section to 0.373 pb with a K factor of about 1.21 for 
the default m c — 1.5 GeV and /i — 2m c . By considering 
its dependence on the charm quark mass and renormal- 
ization scale, the NLO cross section ranges from 0.294 to 
0.409 pb. Furthermore, it will be enhanced by another 
factor of about 1.29 and reaches 0.482 pb when the feed- 
down from ip' is considered. By using da NL (/i) /dfi = 0, 
a maximum a NLO — 0.41 pb is obtained when \i = 1.7 
GeV is chosen, which means the calculation is right in 
the perturbative region. Together with the K factor 
1.21, it could be expected that higher order corrections 
are smaller. By considering theoretical and experimen- 
tal uncertainty, our results are roughly consistent with 
the measurement 0.6 pb from Belle [3, 0], and further- 
more agree well with their recent experimental measure- 
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FIG. 4: Polarization parameter a and angular distribution 
parameter A of J/ip as functions of p with m c = 1.5 GeV and 
/i = 2m c . 



ment 0.43 ± 0.09 ± 0.09 pb [3l|. Thus there is little 
space left for the color-octet contribution now. There 
are large uncertainties in the calculation for the polar- 
ization via color-octet states when the color-octet states 
hadronize into color-singlet states. In contrast, it has 
much lesser uncertainties in the calculation for the po- 
larization of color-singlet state. In the photoproduction 
and hadronproduction of J/ip, there are large discrep- 
ancies between the NLO theoretical predictions [ID, [13] 



and experimental measurements for J/ip polarization. It 
may due to the large contribution from color-octet states 
with large uncertainties, or other mechanism (32J. But in 
this process, color-singlet contribution is dominant and 
the convergence of perturbative QCD expansion works 
very well. Therefore, the prediction for its polarization 
distribution at NLO is well defined and should fit well 
with the experimental data. Our further work[33] gives 
various distributions for e + e~ — > J /ip + X(cc). In order 
to clarify the situation, we suggest to perform further ex- 
perimental analysis of the data based on nowadays huge 
integrated luminosity at the B factories. It is desirable 
that the inclusive J/ip production could be separated into 
e + e~ — > J/ip + X(cc) and e + e~ — > J/ip + AT(non cc) 
in experimental measurement, so that their angular and 
momentum distributions of J/ip polarization can be com- 
pared with their theoretical predictions separately. It 
may be worthwhile to include rclativistic correction ef- 
fects in order to sharpen the test of NRQCD. 

While this paper is being prepared, we are informed 
of the same process also being considered by Ma, Zhang 
and Chao [3J|. Our results are in agreement with theirs. 
In addition, we calculated the momentum distribution 
of J/ip polarization which is a very important issue to 
clarify the J/ip polarization puzzle. 
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